Einstein-Podolsky-Rosen steering nonlocality in pulsed optomechanics and 

macroscopic atomic ensembles 
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A specific scheme for creating Einstein-Podolsky-Rosen (EPR) entanglement between a mechan- 
ical oscillator and an optical pulse has been presented recently by Hofer et al [Phys. Rev. A 84, 
052327 (2011)]. Recent experiments have reported detection of EPR entanglement shared between 
spatially separated macroscopic atomic ensembles [Krauter et al, Phys. Rev. Lett. 107, 080503 
(2011) and Julsgaard et al, Nature 413, 400 (2001)]. We adapt the theories developed for these 
types of experiments, and show how to optimise parameter regimes and measurement strategies, 
so that the EPR steering nonlocality (EPR paradox) can also be realised. For the first type of 
experiment, we show that asymmetric criteria will greatly improve prospects for detecting both 
entanglement and steering. For the second type, steering is predicted to be within experimental 
reach, but modification of the measurement scheme is required to ensure the possibility of spacelike 
separated measurement events. 



I. INTRODUCTION 

The Einstein-Podolsky-Rosen (EPR) paradox [l[ has 
been quite extensively explored for optical fields, both for 
continuous variable (CV) 0-0| and discrete [H, § mea- 
surement outcomes. The evidence for the EPR nonlocal- 
ity is convincing, with some experiments obtaining EPR 
correlations for high detection efficiencies, though not yet 
also for spacelike separated measurement events. The 
EPR paradox has been confirmed for single photon pairs 
[Icj . and mesoscopically, for the quadrature phase am- 
plitudes @ and Stokes polarisation observables of twin 
optical beams [ll| and optical pulses [p3 |. 

The EPR paradox reveals incompatibility between the 
assumption of local realism, and the completeness of 
quantum mechanics. As an example, consistency would 
demand that the physical quantities associated with non- 
commuting observables are predetermined, prior to mea- 
surements, more than can be explained by the Heisenberg 
uncertainty principle. Recently, work by Wiseman, Jones 
and Doherty (l3l . [p4T ] has revealed that the EPR paradox 
is a realisation of quantum steering, a form of nonlocality 
identified by Schrodinger [lfjj where measurements made 
by one observer at a location A can apparently "steer" 
the state of another observer, at location B. Following 
Pil and [l?]], we therefore use the term "EPR steering" 
(or "steering" for simplicity) throughout this paper to re- 
fer to this particular kind of nonlocality, i.e. the EPR 
paradox. 

EPR steering can thus be viewed as a strong form of en- 
tanglement, for which a certain type of nonlocality man- 
ifests. Steering is generally easier to verify than Bell's 
form of nonlocality that falsifies local realism, though, 
in terms of understanding what theories will correctly 
describe nature, may be less informative. Despite that, 
the EPR paradox and steering have been shown to have 
a special usefulness to quantum information tasks (l8r - 
[2l| . including one-sided device independent cryptogra- 
phy [22[ and secure teleportation (23|. For these reasons, 



both fundamental and applied, and perhaps because of 
steadily improving technology, there has been an escala- 
tion in the amount of theoretical interest pjl l24l - |29| and 
in the number of experiments reporting an EPR steering 
nonlocality (Bl. UtI l30l - [32| . In fact, recent experiments re- 
alise EPR steering for photonic qubits, without detection 
efficiency loopholes [33, [34 1, and for spacelike separated 
measurement events [351 ] . 

The evidence however for steering with massive par- 
ticles is almost nonexistent. The Bell nonlocality corre- 
lation strength has been verified without detection effi- 
ciency loopholes for ions [36[ and for qubits in solid state 
systems [371 ] . but not for significant spatial separations. 

We are therefore motivated to analyse the possibility of 
realising steering for systems of massive particles. There 
have already been a number of proposals to detect EPR 
or Bell nonlocality utilising cold atoms or Bose-Einstein 
condensates (BEC) [HI] . These include proposals to de- 
tect C V EPR steering using four wave mixing [3j| , or us- 
ing static and dynamical processes involving atoms con- 
fined to an optical lattice [40l-|42|. Twin atom beams 
[431 . l44j with EPR-type correlate d q uadratures have been 
created in a recent experiment [45j , in a process of pair 
creation that is similar to that used in optics, although 
EPR steering has not been reported. 

In this paper we examine two alternative physical sys- 
tems not yet studied for EPR steering, but that enable 
realisation of quantum limited noise levels. These are op- 
tomechanical oscillators and macroscopic room temper- 
ature atomic ensembles. We analyse the theory of two 
specific experiments, that either demonstrate, or have 
the potential to demonstrate, an EPR-type entanglement 
(EPR-type entanglement is that entanglement existing 
between two spatially separated systems, so that EPR- 
type correlations can be realized), and propose how the 
steering nonlocality might be generated and verified, in 
each case. Both theories model realistically the decoher- 
ence that arises due to thermal reservoirs, and we show 
in both cases that steering is predicted possible over the 
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parameter range considered feasible for the generation of 
entanglement. 

First, in Section HI, we study the recently proposed 
scheme of Hofer et al [46|], to generate EPR-type entan- 
glement between the quadratures of a mechanical oscil- 
lator and an optical pulse, using the radiation pressure 
interaction between the optical cavity mode and the cav- 
ity mirror. A similar proposal was put forward originally 
by Giovannetti et al [43|. We adopt the theoretical model 
developed by Hofer et al, and give predictions for steering 
over the parameter range examined by them. Clearly, the 
realisation of steering is a greater challenge than obtain- 
ing an entanglement. However, for this system in which 
thermal noise plays an asymmetric role, we show that 
the use of asymmetric entanglement and steering crite- 
ria @, 0] will greatly increase the possibility of observ- 
ing both entanglement and steering. One-way steering 
[24 . l3l| is predicted for all values of squeeze parameter, 
while two-way steering becomes possible above a thresh- 
old value. 

Second, in Section IV, we examine the experiments of 
Julsgaard et al [i^ , Krauter et al [5(| , and Muschik et al 
[5l| involving two room temperature macroscopic atomic 
ensembles, for which an EPR-type entanglement has al- 
ready been measured. In this case, modification of the 
experiment is necessary to ensure at least the possibil- 
ity of two spacelike separated measurement events. The 
reported correlation between the macroscopic atomic en- 
sembles is not sufficient to demonstrate steering, but our 
analysis of the theory developed by Muschik et al (Fil [52| 
reveals how to optimise for this form of nonlocality. EPR 
steering is predicted possible within the limits of the pa- 
rameters reported for the experiments, provided one can 
employ local detection of the EPR observables, so that 
conditional measurements can be realised. 



Here Xj are the possible results for a measurement Xj 
at site j where j = A, B, (Xj)\ is the expected value 
of Xj for a given set {A}, and P(A) is the hidden vari- 
able probability distribution function. Nonlocality occurs 
when this LHV model is falsified. Bell nonlocality is the 
case when there is no other constraint placed on the lo- 
cal hidden states. Bell's nonlocality reveals a failure of 
all LHV theories |T]) to describe the measured statistics. 

Subsets of LHV theories can be created in which one, 
or all, of the local hidden states associated with a set of 
spatial locations are constrained to be consistent with the 
predictions of a quantum state. In this case, the model 
CO) is called a local hidden state (LHS) model (iH El E(| . 
We demonstrate a steering and entanglement nonlocality 
when a failure of the LHV is established, where the re- 
striction applies to one, or all, sites respectively [HI, . 
The definitions are qualified by noting that in all cases, 
the nonlocality is to be demonstrated by way of spacelike 
separated local measurements at each of the sites. 

In the bipartite case, the LHV model ((1) applies. The 
model is a quantum separable model, if both local hid- 
den states are constrained to predict quantum statistics, 
meaning that the local averages (Xj)\ (j = A,B) are 
each consistent with a quantum state prediction i.e. can 
be generated by a local quantum density operator. En- 
tanglement as a nonlocality is verified when the quan- 
tum separable model is falsified. Steering occurs in the 
asymmetric case, where the LHS model that is falsified 
assumes only one local hidden state to be constrained to 
predict quantum statistics. It was shown in Refs. (l3l.[T^| 
that the quantum states demonstrating Bell nonlocality 
are a strict subset of the states demonstrating steering, 
which are again a strict subset of all entangled states. 

B. Continuous variable EPR steering criteria 



II. THE EPR STEERING NONLOCALITY: 
DEFINITIONS AND CRITERIA 

A. Definition and the LHS model 

The work of Wiseman et al (HI showed that there are 
three distinct types of quantum nonlocality, depending 
on the restriction placed on the local hidden variable 
(LHV) theory involved to model the measured statistics. 
We consider bipartite nonlocality, where local measure- 
ments are made on two spatially separated subsystems, 
identified as A and B. The observers making measure- 
ments at these two locations are often called Alice and 
Bob respectively. The LHV model introduced by Bell 
[53| identifies two local "hidden" states, one for each site. 
In this model, the predictions for moments of the joint 
measurements made at the two spatially separated sites 
are given by 

(X A X B ) = J^d\P(\)(X A ) x (X B ) x . (1) 



We focus in this paper on CV measurements, where the 
outcomes of measurements are continuous variables. We 
consider two conjugate observables, Xj and Pj at each 
site j = A,B, scaled such that AXjAPj > 1/4 (we are 
using normalised quadratures with |([Ab,Pb])| = 1/2). 
The observables might correspond to a scaled position 
and momentum, or field mode quadrature phase ampli- 
tudes. It has been shown that the Einstein-Podolsky- 
Rosen paradox (48| and steering ('EPR steering") of B 
by A [l| EE Ell is realised if 

A mfA X B A mfA P B < 1/4, (2) 

where A\ nfA X B = A^(X b \O a ) and A\ nf A P B = 
A 2 (P B \0' A ) are the variances of the conditional distri- 
butions P{X B \0 A ) and P{P B \0' A ). Here A) 0' A are 
arbitrary observables for system A, usually selected to 
minimise the variance product [E El Hi- For the EPR 
experiments, the choice is usually either O a = X A and 
0' A = P A , or the reverse choice, depending on the cor- 
relation of the systems. It has been shown that the cri- 
terion @ is necessary and sufficient for steering in the 
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case of the so-called Gaussian states and Gaussian CV 
(quadrature phase amplitude) measurements [l3j (pro- 
vided the observables for Alice are suitably selected). In 
this paper, we examine the predictions for the CV steer- 
ing criterion @, although the LHS model has been used 
to derive other steering nonlocality criteria, including for 
qubits , qudits multi-observable [l7], HH and mul- 
tipartite scenarios |25l . |29| . 

We note that since x 2 + y 2 > 2xy for any real numbers, 
the steering criterion (151) also implies the following sum 
criterion for steering Q. EPR steering of B by A is 
confirmed if 

Ai />A X B + A 2 nLA P B < 1/2 (3) 

since, using the identity, ([3]) implies This sum cri- 
terion, which can also be derived directly from the LHS 
model of [l3j] (see Appendix I) , will be useful to us in Sec- 
tion III, to give predictions for steering in the presence 
of mechanical decoherence. 

We note the inherent asymmetry of the steering cri- 
terion Eq. ©. The asymmetry reflects the asymmetric 
nature of the original EPR paradox, in which it is the re- 
duced noise levels of Alice's predictions for Bob's system 
that are relevant in establishing the paradox [l| . We will 
show that the use of an asymmetric criterion becomes a 
particularly important consideration, even for detecting 
entanglement, where decoherence affects the two systems 
asymmetrically 0, Hl[ • 

C. Symmetric entanglement and steering criteria 

EPR-type entanglement between two modes repre- 
sented by A and B can be verified experimentally us- 
ing the criterion proposed by Duan et al [55j|. For the 
original EPR gedanken experiment [y, a correlation ex- 
ists between Xa and X B , and an anticorrelation exists 
between Pa and P B . In this case, the entanglement is 
verified if the sum and difference variances are reduced, 
according to 

A ent = A 2 (X A + X B ) + A 2 (Pa - P B ) < 1. (4) 

The exact choice of observables for the A field depends 
on the details of the correlation of the EPR system. For 
example, the observation of 

A ent = A 2 (X A + Pb) + A 2 {P A +X B )<1 (5) 

is also sufficient to confirm entanglement between the two 
spatially separated systems. We use the notation A ent 
in either case, to imply an observation of entanglement 
via a measurement of a reduced noise level. It has been 
noted d, [HI that the measurement of 

A ent < 1/2 (6) 

is also sufficient to indicate steering. In fact, because of 
the symmetry in relation to the systems A and B for the 
criteria (H} and ([5]), the observation of A ent < 1/2 will 
indicate two-way steering (that A steers B, and also that 
B steers A) [3l|. 



D. Asymmetric entanglement criteria 

Where there is asymmetry in the system, the sym- 
metric criteria (^flS]) become less sensitive to entangle- 
ment. Other witnesses of entanglement, including nec- 
essary and sufficient measures, have been analysed for 
Gaussian states [5(1 [57}, but we will focus here on the 
local uncertainty relations (LUR) entanglement criteria, 
for which entanglement manifests as a reduction in vari- 
ance caused by an EPR-type correlation. As outlined in 
the Appendix A, entanglement between modes A and B 
is confirmed when 

A 2 (X b -9 x Xa) + A 2 (Pb+ 9p Pa) < i±f^,(7) 
or when 

A 2 (XB+g' x PA) + A 2 (P B +g' p X A ) < l + 9 ^ x9 \ (8) 

where g x , g' x , g p , g' p are arbitrary real constants that can 
be chosen to optimise the entanglement measure, i.e. to 
minimise the variance sum combination. 

Other useful entanglement criteria follow by consider- 
ing variance products [57], HU (see Appendix A) . Entan- 
glement is verified if 

A 2 (XB+g x PA)A 2 (PB+g P X A ) < ( |g * gp J + ■ (9) 

This criterion is useful, since if © holds, the sum crite- 
rion ([5]) will also hold and is equivalent to 

A 2 (X A +h x P B )A 2 (PA+h p X B ) < ( ^M±i y , (io) 

where h x = l/g x and h p = l/g p , which implies a symme- 
try of entanglement detection between the two systems. 

III. EPR STEERING IN PULSED 
OPTOMECHANICS 

There have been proposals to realise EPR correlations 
between a mechanical oscillator and an optical mode j!?} . 
In these proposals, an EPR correlation is generated be- 
tween the position and momentum quadratures X B and 
P B of the oscillator, and the quadrature phase amplitudes 
Xa, Pa defined for the optical field. Recent proposals 
suggest to employ pulsed optomechanics to prepare and 
verify quantum states [59| . and a specific scheme for en- 
tangling a mechanical oscillator and an optical pulse that 
incorporates sources of decoherence has been presented 
by Hofer et al (HWAH) [46(. We examine this scheme, 
in order to understand how it may be modified to detect 
steering. 

A blue detuned light pulse is input to an optomechan- 
ical cavity mode and interacts with the oscillator mirror 
mode via radiation pressure [60| . as described by Hofer 
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et al, and depicted in the Fig. [TJ The boson creation 
and destruction operators for the optical and mechani- 
cal modes are a c , a J and a m , a\ n respectively. Quadra- 
ture phase amplitudes X c i m , P c / m are defined according 
to a c = k(X c + iP c ) and a m = k{X m + iP, n ), where k 
is chosen to ensure the normalised EPR inequality ([2]). 
The key physical parameters are the interaction strength 
<7, the oscillation frequency uj m and dissipation rate 7, 
the optical cavity resonance frequency w c and decay rate 
k, the pulse carrier frequency u>\ and duration time r, 
and no the initial occupation of the thermal state of the 
mirror. 



A. Generation of EPR-type entanglement 

The effective Hamiltonian [HD, [6(| for the system in a 
frame rotating at the laser frequency is 

H = w m o4a m + A c a\a c + g{a m + a\ n )(a c + aj), (11) 




Figure 1: (Color online) Measurement of the EPR entangle- 
ment and steering between an oscillator and a pulse. Step 
I: Entangling a pulse with a mechanical oscillator, (a) Fol- 
lowing HWAH, the blue-detuned pulse interacts with an op- 
tomechanical system. The output pulse quadrature phase 
amplitudes X° ut , P° ut are EPR correlated with the final 
quadratures X™*, P™* °f the mechanical oscillator, accord- 
ing to X° ut Pr* and P° ut X™*, in the limit of a 

large squeezing parameter r. Step II: To verify the steering. 
(b) The output pulse amplitudes X° ut , P° ut are measured 
by homodyne detection. The quadratures of the oscillator 
are measured by interacting the cavity with a second red- 
detuned pulse. The outputs X% ut , Pg ut are a direct measure 
of the original oscillator quadratures, X^ 11 *, P^ 11 , according 
to X% ut = -P™* and Pg uf = X™\ in the appropriate limit 
which optimises the measurement process. To measure steer- 
ing, and to optimise measurement of entanglement, the dif- 
ference currents are combined asymmetrically, by introducing 
an adjustable relative gain factor g. 



where A c = uj c — ui is the detuning of the cavity with 
respect to the laser (4(|. The term in g describes the 
linearised optomechanical coupling due to the radiation 
pressure, and comprises a beam splitter-type coupling 
interaction and a two-mode squeezing type interaction 
[HI, , that is known to generate entanglement [48[ . 

We now summarise the approach of HWAH. The pulse 
shape is assumed to be flat-top and to possess a smooth 
head and tail. HWAH propose the pulse to be either blue- 
detuned, or red-detuned to the cavity resonance (4(|, to 
enable the choice to enhance either the two-mode squeez- 
ing interaction term, for the creation of entanglement, or 
the beam splitter-type interaction term, for the purpose 
of measurement. 

Entanglement is created using a blue detuned pulse. 
In the case where g -C K <C u> m , and after neglecting 
mechanical decoherence effects, HWAH derive a set of 
idealised Langevin equations for the mode operators. To 
justify neglecting mechanical decoherence, they assume 
the pulse duration and time taken for the entire protocol 
to be very short compared to the mechanical decoher- 
ence time. For the blue detuned pulse, after using the 
parameter restrictions to justify a rotating wave approx- 
imation (RWA) and an adiabatic solution for the cavity 
mode, the simplified Langevin equations lead to solutions 
for quadratures X° ut , P° ut of the output pulse, in terms 
of the input pulse quadratures X™, P™ and the initial 
(input) quadratures A™, P™ of the mechanical mode. 
The solutions are 

X° ut = -e r Xl n - Ve 2r - IP™, 

P° ut = -e r P l c n - Ve 2r - IX™, 

X™* = e r X™ + Ve 2r - lP c m , 

C' = e r P™ + Ve 2r ~ IX™, (12) 

where A™* and P™* are the final quadratures of the 
mechanical oscillator, and r = g 2 r/K is the squeezing 
parameter that is dependent on the strength of the inter- 
action and the duration of the pulse. The solutions in the 
limit of large r become X° ut = -e r (X l c n + P™) , P° ut = 
-e r (P™ + X%), A™* = ~P° ut and P™' = ~X° ut . We 
use the subscript c for the pulse quadratures, for conve- 
nience, and stress the distinction between these input/ 
output pulse quadratures and the internal cavity mode 
quadratures (4a. [r33|. 

The input / output solutions allow ready calculation of 
the entanglement and steering criteria of Sec. II. The 
result for the entanglement criterion §5§ is (4f| 

A ent = A 2 (X°f + P° ut ) + A 2 (P^ lt + X° ut ) 

= (n + l)(e r ~Ve 2r -l) 2 , (13) 

where hq is the initial occupation of the oscillator. Our 
choice of scaling means that entanglement is obtained 
when A ent < 1. Further, because we wish to distinguish 
entanglement from the steering, we use the subscript 
"ent" to denote entanglement. In fact, A ent = \Aepr 
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where Arpb is the entanglement measure specified in 
the Ref. j46|. For this system, the EPR correlation is 



between XT* and P°^\ 



and P° ut and X% ut . 



The solutions (flU)) reveal that for large r, A ent — >• 
(no + l)e~ 2r — > 0. This suggests that, potentially, an 
arbitrary amount of EPR correlation, and hence steer- 
ing, can be obtained. The effect of the initial occupation 
number no will be an important factor in determining 
the feasibility of the experiment, however, since there is 
always a practical limitation on the value of r that can 
be achieved experimentally. Figure [2] plots the minimum 
value ro of r required for a given excitation no, revealing 
the logarithmic dependence 



r > rg = In 



n 



1 



2VWTT 



■ In no 



(14) 



discussed in Ref. (46|. We next study the predictions for 
the asymmetric entanglement measures (|81 ll0[) . to show 
an optimisation is possible, so that entanglement can be 
detected for lower values of r, for any given no. 



B. Detecting entanglement 

There is asymmetry in the optomechanical system, 
since the occupation number n cannot be assumed neg- 
ligible for the mechanical oscillator as is justified for the 
optical mode. For this reason, the asymmetric optimized 
criteria (fSl TTU)) become extremely useful for detecting en- 
tanglement in the oscillator-pulse system. 

Entanglement between modes A and B is confirmed by 
(O when 



A 



g,ent — 



The A„ 



A 2 (X^ + g x Pr f ) + A 2 (P° Mt + 9 P X° ut ) 
(1 + g x9p ) /2 



< 1. 

(15) 

±g,ent can be minimized by th e op timal factor 
9x, 9 P via the variational method 1,113, by set- 



tn e op tima 

I El Si, 

ting dDg tSum /dg x . p — 0. For light initially in vacuum 
A 2 Xi n = A 2 P™ = 1/4 and the mirror in a thermal 
state with A 2 X™ = A 2 P™ = n /2 + 1/4, we find the 
optimal choice is 

9x = 9 P 



n /(n + 1) + V /[n /(no + l)] 2 +4e 2 -(e 2 ''-l) 



2e r Ve 2r - 1 



(16) 



Entanglement can be also confirmed by the product 
criterion ([9]) 



A V m 1 ^-i c / V m 1 -J P c > -i 

<+g,prod — 2 ^ L > 

[(\9x9p\ + 1) /4] 

(17) 

which can also be minimized by selecting an optimal fac- 
tor g x , g p . We find the optimal values of g x and g p are 
the same (Eq. (|16p ) as those given for the sum criterion. 
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Figure 2: (Color online) Entanglement detected using the 
symmetric criterion A en t of Eq. (|13[) . Entanglement between 
the oscillator and the optical pulse is detected if A mt < 1. 
Here, A en t = ^Aepr where Aepr is the entanglement mea- 
sure specified by HWAH [46|. (a) A ent is plotted versus the 
squeeze parameter r for various values no = 0, 5, 10, 50, 
where no is the initial mechanical oscillator occupation num- 
ber, (b) The minimum strength ro of the squeezing parameter 
required, for a given no in order to detect entanglement via 
the criterion A ent < 1. 



The entanglement witnesses A 



g,ent, Ag, P rod are plotted 
versus r, for various values of initial occupation number 
no of the oscillator, in Fig. [3l We find the criteria are 
symmetric with respect to interchange of the two sys- 
tems, provided the reciprocals of the values of g x and 
g p are taken, as explained in Section II. The asymmet- 
ric criteria with optimal gain factors, both in the sum 
and product forms, are more sensitive to entanglement 
than the symmetric entanglement witness A ent , which is 
restricted to consider g x = g p = 1. 

We note there is a substantial difference between the 
predictions in the symmetric and asymmetric case. Un- 
like the symmetric criterion (which showed a logarith- 
mic dependence ro = \ In no of the minimum ro required 
for entanglement detection with no, for large no), the 
asymmetric criteria involving A g ent or A g pro( i are not 
sensitive to no, for large no. For a given no, we can al- 
ways show entanglement for any r > 0, provided one can 
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to use a linear estimate <7xP? Mt : based on the result P° ut 
for measurement at A fa [43, li^ . We find 



(pr t ^r t »/2A 2 p c ° 



is 



where = -«X™*,P C ° 

optimised to minimise the variance A 2 n j(X™*|P c tmt ). 
Comparing to the definition of conditional variances 



i\? nf A X B = A 2 (X b |Oa) introduced under Eq. (0, here 
we have selected p° ut as observable Oa-, to obtain the 



conditional variance A 2 n j:(X, 

Similarly, 
evaluated as 



it | pout^ 

Similarly, the conditional variance A 2 



inf 



out\ 



IS 



inf v m 



where g p = X° ut ) + (X° u \ P° ut ))/2A 2 X° ut is 

optimised to minimise the conditional variance. 

The steering of the mechanical oscillator by the pulse 
can be confirmed when (we use the EPR criterion <j2j) ) 



E steer {M\C) =4A, 



if{X ^\P^)A inf (P^\X° c ut ) < 1. 

(20) 

For light initially in vacuum A 2 X" 1 = A 2 P c m = 1/4 and 
the mirror in a thermal state with A 2 X™ = A 2 P™ = 
no/2 + 1/4, we can calculate the prediction for steering. 
We find 



Figure 3: (Color online) Entanglement detected using the 
asymmetric criteria (|7I10|) . (a) A 9jent (Eq. (jTSJ) (b) A g:Pro d 
(Eq. I|17|l) are plotted versus the squeezing parameter r for 
values no = 0, 5, 10, 50 giving the initial occupation number 
of the oscillator. Entanglement between an oscillator and a 
pulse is observed when A g>e nt < 1 or A g4lro d uc t < 1. The 
optimal g x (g p ) to minimize A 9ien t and A 9iPro( 2 is shown in 
the inset of (b). The results indicate improved detection of 
the entanglement, for large no. 



select the optimal choice of gain factor g x , g p . The min- 
imum value ro required for entanglement detection via 
the asymmetric criteria is not (directly) limited by no, 
but will depend in practice on the accuracy achieved for 
selecting the gain factors, which become large, for the 
smaller r values in the high uq limit. 



C. Detecting steering 

We next evaluate the steering predicted by the output 
solutions (fT2|) . The linear relations of the output with 
respect to the inputs ensure Gaussian conditional distri- 
butions [||, and the EPR variances are readily evaluated 
based on an accuracy of inference, that an observer at 
A measuring the output pulse quadratures can predict 
the result for the quadrature X™* of the oscillator, to a 
certain measurable level of uncertainty. A simple way to 
determine this uncertainty for Gaussian distributions is 



A 2 inf (X^\P 



out\ 



via dA 2 nf (X : 



(e r -g x Ve 2r -l) 2 A 2 X™ 

+ (Ve 2 '' - 1 - g x e r ) 2 A 2 Pf -(21) 



which can be minimized by optimizing the factor g x 
™ l \P° ut )/dg x = 0. In fact, we find g x = 

e«- - T(n + 1)] / [e 2r + (e 2r - l)(2n + 1)] . Sim- 
ilarly, the conditional variance A 2 n ^(P°"'|X° Mt ) can be 
minimized by same optimal factor g p = g x . 

Rearranging Eq. (|20|) . we find that the mechanical 



(22) 



2 [e r Ve 



oscillator is steerable by the optical pulse when 



^ T steer 



1 2n 
- In 

2 n - 



1 



1 



1 



In 2. 



For this case, the requirement on the minimum strength 
of effective optomechanical interaction (for a steering of 
the mechanical oscillator by the optical pulse) approaches 
a constant, r ste er = ^ In 2 = 0.35, as the initial occupa- 
tion no of the oscillator increases, i.e. as no — > oo. Figure 

plots the values of E s teer(M\C) given by Eq. ([20)) ver- 
sus r, for various values of occupation no- 

The steering of the pulse by the mechanical oscillator 
can be certified if 



E steer (C\M) = 4A m/ (xr t |Pr t )A m /(P ( 



out | J£OUt\ 



') < 1. 
(23) 

We find that the optical pulse is steerable by the mechan- 
ical oscillator for any n and r > 0. Figure 2Jd presents 
the results for different no- For a squeezing parameter 
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Figure 4: (Color online) Detecting steering between the oscil- 
lator and the pulse, (a) The steering nonlocality parameter 
E a teer{M\C) (Eq. (J20j) ) of the oscillator M by the cavity C, 
is plotted versus r for various values of the initial oscillator 
occupation no. Steering is detected when E s teer(M\C) < 1. 
Strong steering implies E s teer — > 0. The inset shows the opti- 
mal g = g x = g P to minimize E a teer(M\C). (b) Certification 
of one-way and two-way steering versus the squeezing parame- 
ter r. The conditional variance products according to criteria 
(12011 and (|23[) are shown for various no. 



r < ^ln2, one-way steering is predicted, whereas for 
r > \ In 2, two-way steering becomes possible, mean- 
ing that both mirror and pulse can steer the other re- 
mote subsystem. Clearly, the asymmetry of the steering 
is due to the asymmetry of the thermal effects on the 
two systems. As might be expected, the presence of the 
thermal mechanical excitation no makes the steering of 
the oscillator, by the pulse, a more difficult challenge 
experimentally. This result is however of greater funda- 
mental significance, in view of the original issues raised 
by Schrodinger [l3j], since in this case, it is the massive 
mesoscopic system that is being "steered", by the distant 
observer (the pulse system). 

We find a similar prediction is given for paramet- 
ric amplification in the presence of asymmetric ther- 
mal noise, as modeled by the interaction Hamiltonian 
Hj = ihK(a'w — ab), for two modes with boson opera- 
tors a and b. Generally, the presence of additional ther- 



mal noise in mode b creates the possibility of one-way 
steering, making steering of 6 by a more difficult than 
steering of a by b. This contrasts with the result for dis- 
sipation. When mode b has additional losses, steering of 
a by b is more difficult than b by a @, 

An application of EPR-type entanglement is to enable 
teleportation of a quantum state [64| [65[ . HWAH present 
a teleportation protocol based on the realisation of con- 
tinuous variable teleportation using optical fields. The 
teleportation fidelity is given by the overlap of the re- 
constructed state p ou t and an arbitrary quantum state 
\tpin)- F = (tpin\pout\ipin) [H- For coherent input states, 
and where systems are symmetric, the fidelity is given 
by the entanglement parameter of Duan et al. [55j . 
F = 1/(1 + A ent ). The value of F > 1/2 (A ent < 1) 
is required to do better than an optimal classical strat- 
egy [60|. However, to teleport with security, so that the 
quality of the state transported is sufficiently high that 
there can be no replica states of the same or better qual- 
ity at another location, higher fidelities are necessary [23| , 
F > 2/3. This requires A ent < 1/2 in the CV case, the 
level of correlation ^ given by symmetric steering. 



D. Mechanical decoherence 

The work of Ref. (4(| also studies the predictions for 
A ent in the presence of mechanical decoherence. Their 
work presents predictions showing A ent < 0.5, which is 
sufficient to demonstrate the steering (by ©), for large 
thermal bath excitation numbers n, provided the me- 
chanical quality factor Q m = Lu m /"f is high, of order 
~ 10 5 , and provided the oscillator is precooled to the 
ground state (no ~ 0). In this case, the detection of the 
steering is clearly more challenging, though in principle 
feasible. 

The results from the previous section however indicate 
a more promising prediction if the asymmetric entangle- 
ment and steering criteria are used. Here, the sensitivity 
of the entanglement to the initial occupation no of the os- 
cillator is greatly reduced, and detection of entanglement 
and steering becomes possible at higher no, for lower, 
more accessible values of squeeze parameter r. This re- 
sult changes the implication for the effect of mechanical 
decoherence by broadening the available parameter range 
available for detection of EPR-type entanglement. 

To give an indication of this, we have analysed the 
detailed model presented by HWAH, for the feasibility 
of detecting entanglement and steering using the asym- 
metric criteria. Details are given in Appendix B. The 
model focuses on the regime of parameters for which 
j«k<C Lj m . which would justify the idealised approach 
given by them that is described in the previous Sec- 
tion. HWAH define dimensionless parameters rj = K/uj m , 
£ = g/n, and e = 7T where 7 is the decoherence decay 
rate, and use their model to arrive at optimal choices of 
these parameters, for the detection of A e „ t in the pres- 
ence of the mechanical heat bath with occupation number 
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Figure 5: (Color online) Entanglement and steering in the 
presence of mechanical decoherence. (a) The entanglement 
Aent and A a . ent (Eq. I24[) including the effect of mechanical 
decoherence. Here n is the thermal occupation of the associ- 
ated bath. The A 9jen t is plotted with the optimal gain factor 
(fl — 9x = 9p), whereas A en t corresponds to the choice 3 = 1. 
The A g , ent predicts a lower variance for no = 50 (black solid) 
in comparison with A en t (red solid). For no = 0, the optimal 
value is g — 1, and the results for A ent and A St ent agree (blue 
dash-dotted) and (light blue solid). The values of (n + ^)e 
are plotted as the light blue dashed line (no = 0) and the red 
dashed line (no = 50). (b) Detection of two-way steering in 
the presence of mechanical decoherence for the case no = 50. 
Here we plot E„t e er,sum and use the sum steering criteria of 
Eqs. jB5j and (|BTu) . 



n, for an initial excitation tiq. The value of n depends on 
the temperature of the bath. 

HWAH present results using the symmetric sum crite- 
rion A en t. To enable a more direct comparison, we also 
consider the sum criteria, but generalised to allow for 
asymmetry i.e. we consider the asymmetric entanglement 
and steering in the sum format, given by criteria (|15p and 
Q . The entanglement witness (fT5|) is normalised, and is 
written 



A Sient = A S]ent | 7=0 + (2n + l)e/ (1 + g x g p ) , 



(24) 



where A ffjerlt | 7= o is the value with no decoherence. The 
asymmetric entanglement condition A„ j£ , nt becomes A ent 



with the choice g = g x = g p = 1, and this prediction for 
the parameters chosen by HWAH is replicated in Fig. [S£i, 
for the two values of no- By optimizing g, the A g ^ nt can 
be reduced, to allow better predictions for entanglement. 

The values A ent | 7= o = A g , e „ t | 7=0 where g = g x = 
g p = 1 reported by HWAH correspond to those with no 
mechanical decoherence, and in the analysis of HWAH 
are obtained from the full Langevin model. However, 
assuming the parameter range is consistent with g <C 
k <C u m , the solutions A erat L = n and A 



*g,ent I7— 



= n cannot 



deviate significantly from those obtained in the idealised 
model given by the solutions (fT2"|) . Hence, we obtain an 
approximate prediction for the asymmetric criteria by us- 
ing the idealised model to solve for r, given the values for 
A e „t (for the given no). This enables the reconstruction 
of the asymmetric prediction for A ffjen t| 7= o, which for 
the higher value of n$ = 50, is considerably lower than 
A ent | 7= o. We then give a possible prediction for the en- 
tanglement for values considered reasonable by HAWH, 
by selecting the choice e = e op t given in their Fig. 2 
(Appendix B). 

The resulting predictions for entanglement detected by 
the symmetric and asymmetric witnesses is plotted in 
Fig. [Sfe,, and shows the enhanced possibility to detect 
entanglement in the presence of thermal mechanical de- 
coherence, without the need to use laser cooling to reduce 
the value of no, if one uses the asymmetric criteria. The 
analysis for steering is also given in the Appendix B. The 
observation of A ent < 0.5 is sufficient to verify steering, 
and this prediction can be read directly from the results 
already presented in Ref. (4(|. However, we present the 
actual prediction for the steering parameter, optimised 
for g in Fig. [5Jd. As for entanglement, the results indicate 
the potential to detect two-way steering in the presence 
of thermal mechanical decoherence, without the need to 
use laser cooling to reduce the value of no- 



E. Measurement of EPR steering 

Our calculations suggest that steering between a me- 
chanical oscillator and a pulse is distinctly feasible. The 
method of detecting steering must be considered care- 
fully, to ensure a demonstration of quantum nonlocality 

m ■ 

To verify the steering, one needs to measure the 
quadrature phase amplitudes X™* = X° ut and P™' = 
pout Q £ blue-detuned output pulse, at A, and also (in 
principle simultaneously) at a different location, those of 
the final state of the mechanical oscillator, at B. The pro- 
posed scheme of HWAH enables such measurements to be 
made. The output pulse measurements are made with a 
homodyne detector, as depicted in Fig. [TJ A proposal for 
measuring the final quadratures of the mechanical oscil- 
lator is outlined in HWAH. A second red-detuned pulse 
is coupled to the oscillator, so that the output pulse and 
final mechanical quadratures, I™', P° ut and X™\ P 7 ° n ut 



respectively, are given by [46 
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B J 



(25) 



where r' is the interaction strength for the second pulse. 
In the limit of large r' — > oo, the output pulse quadra- 
tures are a measure of the original quadratures of the 
mechanical oscillator: X% ut = -P™ and Pg it = X™. 
With this assumption, measurement of the quadratures 

X out^ pout and X out^ pout of the twQ Qutput pulses A 

and B will enable demonstration of the steering nonlo- 
cality, in accordance with criterion The quadrature 
measurements are performed via homodyne detection, as 
illustrated in Fig. [TJd. 

Since the distributions are predicted to be Gaussian, 
by the theory of Refs. (4(| |U [6(j, the procedure used in 
the experiment Q involving a variable control g would 
enable measurement of the variances of the conditional 
distributions, with only a minor adjustment of the pro- 
cedure for measuring A ent = ^A epr . The method out- 
lined in Ref. [4(| using optical homodyne tomography 
would enable verification of the Gaussian prediction, and 
evaluation of the conditional distributions, for a rigorous 
demonstration of the EPR steering nonlocality. 

For rigorous EPR tests, the measurement time To 
of the mechanical oscillator quadrature amplitudes also 
needs to be considered. This time is limited by the second 
pulse duration and cavity interaction time, which ideally 
would be short enough to enable a spacelike interval be- 
tween the measurement events at A and B. On the other 
hand, the time must be long enough that the interaction 
strength r' enables an accurate measurement of the me- 
chanical quadratures. The use of squeezed pulses would 
enhance the sensitivity of the measurements, by reduc- 
ing the noise in X l £ (Pg 1 ), for the measurements of X out 
and P out respectively. 
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Figure 6: (Color online) Entangling two oscillators: (a) Gen- 
eration of the entanglement takes place when the output of the 
first cavity is injected into the second cavity, as red-detuned. 
The final states of the two oscillators at A and B will be entan- 
gled, (b) The entanglement can be verified, at a later stage, 
using two red detuned pulses, and the homodyne scheme set 
up as depicted, to measure the conditional inference variances 
A inf X mj and A inf Pmj • 



ond cavity. Using relation Eq. (fT2"|) 



j^out 

pout 

c 

Xout 
ml 



-e r X™ - x/e^-LFd, 



2r 



1x1 



e r xz + v^Tp: 



(26) 



as above, we see the first blue-detuned output pulse is 
entangled with the final quadratures of the mechanical 
oscillator (1). In the limit of large squeezing parameter r, 
the relations are = -P° ut and P£f = ~X° ut . The 
output (of the first) pulse is then red-detuned relative 
to and incident on a second mechanical oscillator system 
(m2), so that the input-output relations (|2"5")) hold 



F. Entanglement and steering between two 
mechanical oscillators 

It would be a challenge to demonstrate steering be- 
tween two massive systems. Bipartite steering entangle- 
ment between two mechanical oscillators can be achieved, 
in principle, by swapping the entanglement between the 
oscillator ml at location A and the output pulse at lo- 
cation B, to an entanglement between the oscillator ml 
at A and a second mechanical oscillator m2 at location 
B (Fig. [5]). Thus, we note from the input/ output re- 
lations that one can EPR entangle two oscillators, if the 
first entangling pulse subsequently interacts with a sec- 



X° B ut = ~e- r 'x%-Vl-e-^'P™ 2 , 

P° B ut = -e- r 'P B n + Vl- e- 2r 'X™ 2 (27) 

vout _ -r' yin , f\ „-2r' pin 
A m2 — e A m2 + V t ~ e fg , 

pout _ -r' yin f, _ 2r' V in 

r m2 — e A m2 V 1 c ^B ■ 

In the limit of large r', we find X%$ = Pp and = 
—Xg l so that the entanglement has been swapped into 
that of the mechanical oscillator. The relations be- 
tween the quadratures of the mechanical oscillators is 
*Sa = -X°Z and P^ = P°f as a§ = a?*. We note 
that similar bipartite entanglement could potentially be 
realised by schemes that use two-mode entangled light 
fields (67|. 
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IV. EPR STEERING BETWEEN ATOMIC 
ENSEMBLES 

We next examine the possibility for steering between 
two macroscopic room temperature atomic ensembles. 
For this system, an EPR-type entanglement has already 
been realised using the criterion (U). In this case, there 
is symmetry between the two ensembles, and the asym- 
metric criteria for entanglement prove not to be so use- 
ful. However, to enhance detection of steering, it be- 
comes essential to measure the conditional EPR paradox 
(steering) criterion Eq. @, and to devise a measure- 
ment scheme which enables separate local measurement 
on each system. 



A. Generation of EPR entanglement 

The experiments of Julsgaard, Krauter, Muschik and 
co-workers [49M5H achieve entanglement of two macro- 
scopic spatially separated atomic ensembles. A schematic 
diagram of their experiments is given in Fig. [7^,. Entan- 
glement is detected by measurement of A ent < 1, using 
the criteria of Eqs. Q and (fT5|) [55| . The experiment 
does not yet achieve steering, for two reasons. First, it 
has not been shown that an EPR paradox or steering cor- 
relation criterion of the type discussed in Refs. (39 . HH 
has been satisfied. The noise reduction of A ent < 0.5 
(Eq. ([5])) would be sufficient to demonstrate this, as 
would conditional variance measurements that satisfy Eq. 
©. Second, the EPR observables, which in this case are 
the collective atomic operators, are not each measured 
locally. 

Analysis of the theoretical models given in Refs. 
(49l . loll. |68| suggest that a demonstration of steering is 
within experimental reach. We summarise the theory 
and the associated experiments. The two spatially sepa- 
rated atomic ensembles are denoted by A and B, and 
Schwinger collective spins, J X / B , J^/b- ^a/b are 
fined for each atomic ensemble, which each contain N 
atoms. The operators are defined with respect to two 
atomic levels, denoted, for the jth atom in ensemble A, 
as |l)j and \2)j. Thus we define: jf = ~ 

|2)(2| i )/2; Jf = E7=i(|l)(2|i + |2)(l| J -)A and Jj = 

Y^jLi — |2)(l|j) /2i. Similar operators are de- 

fined for two selected levels of the ensemble B. Each 
atomic ensemble is prepared, using a detuned pump laser 
pulse, in an atomic spin coherent state with a large mean 
spin J x , so that (J x ) ~ N, and this implies the atoms 
are prepared in a superposition of states |1) and |2). The 
mean spins are equal and opposite for the ensembles A 
and B, i.e. (J x ) = ~{Jb)- Based on the Heisenberg un- 
certainty relation, the steering criterion of the type @ 
becomes Q 

A m/ Jf A m/ j£ < I|(J* )|. (28) 




Figure 7: Schematic diagram to test the EPR steering non- 
locality between atomic ensembles, (a) Depiction of the ar- 
rangement of the experiments of Refs. [49h51||. which measure 
the sum collective spins Jjf + Jf or ,l\ + J Y . (b) Depic- 
tion of an experiment to demonstrate the steering between 
two atomic ensembles. Independent local measurements are 
made on the atomic ensembles A and B, by measurement of 
the Stokes observable S Y of the transmitted verifying pulses. 
The choice of whether to measure J z or J Y at each ensemble 
(which is made by a rotation of the atomic spin of the ensem- 
ble) is made independently and after preparation of the EPR 
state (which is done by the entangling pulses). The measure- 
ment process becomes more sensitive if squeezed input pulses 
are used, as described in the text. 



The two atomic ensembles are entangled using a de- 
tuned polarised laser pulse. The laser field is described 
by the Stokes operators S x , S , S , where S x is the 
difference between photon numbers in the orthogonal X 
and Y linear polarisation directions; S Y is the difference 
in number for polarisation modes rotated by ±7r/4; and 
S z is the difference between the photon numbers in the 
two circular polarised modes defined relative to the prop- 
agation direction Z. The Stokes operators are written 
in terms of 011,2, the operators for the circular polarised 
modes. Thus, 




S x = \a\ai + ai4j A 

S Y = (a\ai - a 2 a\) /2i. (29) 

In the experiment, the light is polarised along the x di- 
rection, so that (S x ) ~ N p where N p is the number of 
photons in the optical pulse. 

Following the theory of Duan et al. (68[, the detuned 
entangling pulse propagates through an atomic ensemble, 
and the outputs are given in terms of the inputs according 
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(30) 



The a and /3 are constants, and the subscripts out and 
in denote, for the field, the outputs and inputs to the en- 
semble, and, for the ensemble, the initial and final states. 
The /? is reversed in sign for the second ensemble, so that 
after successive interaction with both ensembles, the final 
output is given in terms of the first input as 
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(31) 



Here the subscripts A, B denote the spin operators for 
the A, B ensemble. As described in the Refs. (iil. l68|. 
provided a is large enough, the field output S Y ut gives a 
measure of the collective spin Jjf + J B , which is a con- 
stant of the motion i.e. not affected by the interaction 
process. The Stokes parameter S Y ut of the output field 
is measured, by polarising beam splitters and detectors, 
and the ensembles prepared in the state for which ideally 
the value of J z + J B is known and constant. In practice, 
the preparation implies a state with reduced noise level 



in J z + J§. The collective spin Jj - 



J B is also a constant 



/t TJ B; LU11CLL1VC OjJlll 'J ^ 

of the motion, and hence the ensembles can be prepared 
via a second pulse (and by rotating the atomic spin) in a 
state with reduced fluctuation in j\ + J B . The reduced 
fluctuation in the noise levels according to 



A,. 



AVi + Jf) + A 2 (Jl + j£) 



< 1 (32) 



signifies entanglement between the two atomic ensembles 
[55| . Details of how the A ent and the mean spin (J A / B ) 
are measured are given in the Refs. j49l . [68j . 

The key point, from the perspective of this paper, 
is that the field-atom solutions Eq. (|3"Tj) give no lower 
bound to the amount of noise reduction in jf + and 



J A 



J Y that is possible. 



A ' B 

B uncLu id jjuoniuic. The solutions Eq. (|30[) on 
which the prediction is based are valid in the limit where 
damping effects can be neglected. Duan et al. [68[ have 
analysed full atomic solutions, and shown this regime is 
achievable when N ~ N p and provided the field detun- 
ings are much greater than spontaneous emission rates. 
Hence, generation of the EPR paradox between the two 
ensembles is predicted possible in principle using this 
technique, since enough EPR correlation will imply steer- 
ing, by ©. 

The experiments |49l . [5fJ detect entanglement by mea- 
surement on a verifying pulse that is transmitted through 
both ensembles. Modification is needed for an EPR ex- 
periment. One possible strategy is to use two verifying 



pulses, defined with Stokes parameters S A 
respectively, one propagating through each ensemble A 
and B respectively (as shown in Fig. [7}d). The outputs 
in terms of the inputs are given by the solutions Eq. pop , 
so that 



qY . qY 



in,Ai 



Sout,B — Sin,B + a Jin,B- 



(33) 



Measurement of S Y ut A at one location, and S Y ut B at the 
other location enables local determination of Jjf and J§, 
as required for a test of EPR nonlocality. The local mea- 
surements of j\ and J B can be made similarly. The 
sensitivity of the measurement of the collective atomic 
spins J A j B would be improved if the input fields are 
"squeezed", so that, for example, when measuring jf and 



J§, squeezed fields with A 5, 
respectively, are used. 



Y 

in, A 



and AS£ B 



0. 



Detailed calculation of steering for an 
engineered dissipative system 



The recent experiment of Krauter et al. |50| employs 
an engineered dissipative process to generate a long-lived 
entanglement between the atomic ensembles. Entangle- 
ment values of A ent ~ 0.9 have been realised experimen- 
tally. We therefore analyse in more depth the theory 
presented by Muschik, Polzik and Cirac (MPC) [H [52j 
for the experiment, to calculate the degree of EPR steer- 
ing predicted for the schematic set-up of Fig. [7)d, given 
the parameters of their experiment. 



The MPC model introduces operators A = p,J, +vJ 



B • 



B 



vJt, where J 



operators with J 



A/B 
1 



1 v^W 
2^7 



denote collective spin 

I t)U \i and j + = 



I |)(t |j such that J Y = (J+ + J~)/2 and 

The p and 



JN ^ = 1 

J z = i(J+ - J-) /2, and p 2 - v 2 = 1. 
v characterise a squeezed state with squeeze parame- 
ter r, where p = cosh (r) and v = sinh(r). The 
correlations are described by a master equation dtp = 
x (Apti - A*Ap + BpB^ - B^Bp + H.c.) + C noise p, 
where p is the atomic density operator, d is the optical 
depth of an ensemble, T is the single atom radiative de- 
cay, and C no i se p describes undesired processes such as 
single atom spontaneous emission noise. The Linblad 
terms given in the parentheses drive the system into the 
EPR state [5(j- Other parameters are the number of 
atoms TV-j- and JVj, in each of the two levels; and the pop- 
ulation P 2 (t) = (N t - N i )/N 2 (t) where N 2 = N t + N± 
is the number of atoms in the two-level system. In Ap- 
pendix C of their paper (Hlj . equations are derived for 
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the evolution of the atomic spins: 

d t ({j z f) = -[f + dTP 2 (t)]((j z ) 2 ) 

+ ^[f + dTP 2 (t) 2 ( P 2 + v 2 )]. 
d t ((j Y ) 2 ) = -[f + dTP 2 (t)](j Y ) 2 ) 

+ ^[f + dTP 2 (t) 2 ^ 2 + , 2 )]. 
dt(JlJ Z B ) = -[f + dTP 2 (t)}(J z J z ) 

+ E fliydr p 2 (t) 2 , 
dt{J Y A Jl) = -[T + drP 2 (t)](J Y J Y ) 

-^ v dTP 2 (t) 2 , 



(34) 



where f = T coo i + T heat + T d , T coo i(T heat ) is the to- 
tal single-particle cooling (heating) rate and Yd is the 
total dephasing rate. This leads to the equation (CI) 
of their paper, describing the evolution of the variances 

To obtain the optimal prediction for steering in the 
Gaussian system, we modify the analysis to calculate the 
evolution of the conditional variances, Var(J A \J B ) = 
A 2 (J Y ±g Y J Y B ) and Var(J z \J z ) = A 2 (J Z ± g z J z ). 
The conditional variances can be measured by the ar- 
rangement of Fig. [7J}, which allow a relative weighting of 
the measurement of the collective spin of one ensemble, 
against that of the second ensemble. The equations are 



d t Var{Jl\J Y ] 



d t Var{Jl\J z ) 



= -[f + dTP 2 (t)}Var{j\\J Y ) 



+ - 



■ [r + drp 2 (t)V + v 2 )] 



TNg Y ^dTP 2 (t) 2 , 

-[T + dTP 2 (t)] Var(J z A \3 z ) 

+ N{1+ A 9z) [T + dTP 2 (t) 2 ^ 2 



±Ng z ^dTP 2 {tf. 
The steady state solutions are given by 



(35) 



Var{Jl\J Y 



b y°o 



A Y [f + dTP| >00 (M 2 + v 2 )] T Ngy^dTP 



2 ,oc 



[f + dTP 2)0 o] 
VariJilJi)^ = 

A z [f + dTP 2 ^ 2 + v 2 )] ± Ng Y ^dTP 2 c 



(36) 



where Ay = N(l + g Y )/^ and A z = N(l + g 2 z )/4:, and 
we will choose 



9y = -gz 



±/iudrP|, 



[f + dTP 2 00 (^+^)]/2 



(37) 
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Figure 8: (Color online) Steady state EPR steering using 
the engineered dissipatively driven system of Muschik et al. 
[5l| . (a) The predictions for the steady state entanglement 
Aent = £oo versus Z = (/^ — «) _1 for an optical depth <f = 30 
per ensemble. We use the parameters presented in Fig. 3 
of Ref. [5lfl . < 1 represents entanglement between two 
ensembles, and the observation of A ent = < 0-5 i s a l so 
sufficient to verify steering, (b) The predictions for the steer- 
ing (|39[) using the same parameters but optimal values of gain 
factor gy and gz- E ste er,oo < 1 indicates two-way steering 
between ensembles. The inset shows the optimal gain factor 
9Y(-gz) for minimizing E at eer,oa- In this case r cooi = fi 2 T 
and T heat = v 2 T. The dephasing rate F d = V r d ad + T a d dd con- 
sists of a radiative part V d ad — 2(/i 2 + v 2 )F which is due to 
light-induced transitions, and an additional term Y d dd which 
summarizes all nonradiative sources of dephasing. The violet 
line is for pure radiative damping, T d dd = 0. T is the single 
particle decay rate. 



to minimize the conditional variances. The steady state 
solution for \(J A B }\ is given by Muschik et al. [51|, as 



\{Ja,b) 



N 



(38) 



for t — > oo. 

The steady state steering E stee r, oo is then given as 



yJVariJ^UVariJflJiX 



(39) 
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The entanglement is 

a _ Var(JX\J%)oo + Var(J%\J§U 

(K J f)L + \9Y9z\\{J§)\J /2 

where we also use introduce the notation ^ used by 
MPC: is the value of A erlt in the steady state. In 
fact, due to the symmetry of the two systems, A g ^ ent is 
minimized by the optimal gain factor Qy = —gz = lj 
and hence agrees with Eq. (15) in Ref. [51j . The steady 
state entanglement ^ = A ent is plotted in Fig. [S^,. 

We plot the optimal steering prediction, for the param- 
eters chosen by Muschik et al. J5l| in their Fig. 3. Of 
course, the observation of A ent = ^ < 0.5 (by ([6])) is 
also sufficient to verify steering, but the noise reduction 
in their case is not enough to establish this. We there- 
fore present the actual prediction for the steering, opti- 
mised for gy{gz), in Fig. [SJd. Steering is clearly more 
difficult to generate than entanglement, with values of 
E s teer.oo ~ 0.9 being predicted for mainly pure radiative 
damping, for which values of entanglement A ent ~ 0.6 
are predicted. Given that the experiment has achieved 
steady state entanglement, a measurement of steering 
would not seem an impossibility, however, even for room 
temperature atoms. 



V. SUMMARY 

In summary, the EPR paradox entanglement studied in 
this paper confirms an inconsistency of local realism with 
the completeness of quantum mechanics, and is evidence 
for the form of nonlocality called "EPR steering". We 
have analysed the possibility of realising such steering 
for two physical systems: pulsed cavity optomechanics, 
and spatially separated macroscopic atomic ensembles. 
We propose how steering might be generated, verified, 
and optimised in each case. 

Especially, we note that for the first system, because 
thermal noise is more significant for the mechanical os- 
cillator than for the optical mode, the use of asymmetric 
entanglement and steering criteria will greatly increase 
the possibility of observing both entanglement and steer- 
ing. This is true even where there is a significant initial 
occupation number of the oscillator and in the presence 
of mechanical decoherence. We find a regime of one- 
way steering is possible, where the optical system can be 
steered by, but cannot steer, the mechanical oscillator. 

For the second system, EPR steering is predicted 
within the parameter range of the experiment, provided 
nonradiative sources of atomic dephasing are not too 
large. It becomes essential however to modify the ex- 
perimental design to ensure detection of two local EPR 
observables, so that conditional measurements can be re- 
alised. 

In conclusion, we have considered detailed theories for 
entanglement generation in two sorts of massive systems. 
The theories model the major sources of decoherence, 



(dissipation and thermal reservoirs) and for the case of 
gas ensembles have been partially verified experimentally. 
Our work reveals that detection of steering is more chal- 
lenging than detection of entanglement, but is predicted 
feasible, for accessible parameter regimes. 
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Appendix A: Proofs 

We outline the derivation of the EPR entanglement 
criteria used in this paper. Assuming separability, and 
that the variance of a mixture must not be less than the 
average of the variances of its components, we can write 

M 

A 2 (A B - g x X A ) + A 2 (P B + g p P A ) 
> Pr(&r(Xb - 9 X X A ) + A 2 R (P B + g p P A )) 

R 

= Y, P R[{X 2 B )R + gl{X 2 A ) R - 2g x (X B ) R (X A ) R ] 

R 

+ E P r[( P bU + 9l(Pl)R + 2g p (P B ) R (P A ) R ] 

R 

- E Pr(X b - g x X A ) 2 R - E Pr(Pb + g P P A ) 2 R 

R R 

= Y, Pb(& 2 rXb + A 2 R P B + g 2 x A 2 R X A + g 2 p A 2 R P A ). 

R. 

(Al) 

Here, the subscript R denotes the variance or average for 
the state depicted by R (namely p A or p R ). Separabil- 
ity implies the factorisation, for example, (X A X B ) R = 
(X A } R (X B ) R . For any quantum state, uncertainty re- 
lations follow from AX B AP B > \([X B , P B })\/2 = 1/4. 
The local uncertainty relation is A 2 A^ + A 2 P B > 1/2 
and A 2 (g x X A ) + A 2 (g y P A ) > g x g p /2 (we use the identity 
x 2 + y > 2xy, for x, y real numbers). The inequality 
thus becomes 

A 2 (X B -g x X A ) + A 2 {P B +g p P A ) > (l+g x9p )/2, (A2) 

which if violated confirms entanglement between the 
modes denoted by A and B. For g x = g p = 1, the crite- 
rion reduces to the symmetric criterion for entanglement 
used by HWAH. 
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We now note that EPR steering criteria can be derived, 
from the above proof, by using the LHS model of Ref. 
[l3j . In this case, the second site that is being steered 
remains with a quantum restriction, meaning that the 
uncertainty relation holds for it, but the steering site be- 
comes unconstrained (except as a local hidden variable 
site). In that case, the right side bound is replaced by ei- 
ther 1/2 or g x g p /2 depending of whether we wish to show 
A steer B, or B steers A respectively. For g x = g p = 1, 
we obtain the criterion © for EPR steering. 

Criteria involving the product of variances can also be 
derived (as in Refs. (57l |58|): 

A 2 {X B - g x X A )A 2 (P B + g p P A ) 

> [J2 Pr&KXb - g x X A )] [J2 Pr& 2 r(Pb + 9 p Pa)] 

R R 

= J2 PrI^rXb + g 2 x A 2 R x A ] x J2 Pr[^ 2 rPb + g 2 P A 2 R P A ] 

R. R 

> {J2 Pr[A%X b A 2 r P b + g 2 x g 2 p A 2 R X A A 2 R P A 

R 

+g 2 x A 2 R X A A 2 R P B + g 2 p A 2 R P A A 2 R X B ]^ 2 } 2 

> ^ + 9 2 X 9 2 P + 2\g x g P \ _ ^ \g x g P \ + l ^ 2 (A3) 

Here we have used the standard Cauchy Schwarz inequal- 
ity and the Heisenberg uncertainty relations. Hence, en- 
tanglement is verified when 

A 2 (X B -g x X A )A 2 (P B +g p P A )< nM^±i y , (A4) 

Appendix B: Mechanical decoherence 

We analyse the model for mechanical deocherence pre- 
sented by HWAH [4(| , for the feasibility of detecting en- 
tanglement and steering using the asymmetric criteria. 
The model focuses on the regime of parameters for which 
g <C k <C uj m . Dimensionless parameters rj = K,/ui m , 
£ = g/n and e = 7T are defined, where 7 is the decoher- 
ence decay rate, and HWAH use their model to arrive at 
optimal choices of these parameters, for the detection of 
A ent in the presence of the mechanical heat bath with 
occupation number n, for an initial excitation no. The 
value of n depends on the temperature of the bath. 

HWAH arrive at solutions (Eq. (21) of their paper) for 
the quadrature phase amplitudes, 

\rout -yout I 1 / jji 

A m — A m | 7 =0 + Ve^s, 

C = Pr°n U %=o + V~eFc (Bl) 

where the additional noise operators F s and F c are due to 
thermal mechanical noise. Thermal noise is neglected for 
the optical fields, for which X° ut - A c °"*| 7=0 - Also, the 
effect of damping is justified to be small compared to the 
effect of the thermal noise term. The model leads, under 



the approximations that justify a perturbative expansion 
for the symmetric entanglement 

A ent = A ent | 7=0 + (n + -)e. (B2) 

We now analyse the prediction for the asymmetrical 
entanglement criterion. We find, using (|B1[) . that 

A 2 (A ? r + g x P° ut ) = A 2 (xr + g x P° ut %=* + \{n + 
A 2 (Pr* + 9 P X° ut ) = A 2 (P r r + 9 P X° U % =0 + |(n + 

( 

because the extra noise arises solely from the mechanical 
quadratures. 

For a direct comparison, we consider the entanglement 
and steering in the sum format, given by criteria (|15[) 
and ([3]). The entanglement witness (fTS)) is normalised, 
and becomes (using (|B3|) ) 

Ag tent = Ag ient | 7=0 + (2rT+ l)e/ (1 + g x g p ) . 

(B4) 

The prediction for the parameters chosen by HWAH in 
their Fig. 2 is replicated in the plot of Fig. [5^ for 
g = 9x = 9p = 1) for the two values of no considered by 
them. By optimizing g, the A g ^ ent can be reduced. The 
values A e „t| 7= o reported by HWAH correspond to those 
with no mechanical decoherence, and in the analysis of 
HWAH are obtained from the full Langevin model. How- 
ever, since the parameter range considered by HAWH is 
consistent with j<Ck« uj m , the solutions A ellt | 7=0 and 
A 9i e„t| 7= o cannot deviate significantly from those ob- 
tained in the idealised model given by the solutions (TT3]) . 
Hence, we approximate the prediction for the asymmetric 
criteria by solving the values given for A en t = Aepr/% 
(for the no) for r. This enables the reconstruction of 
the asymmetric prediction for A gjent | 7= o, which for the 
higher value of no = 50, is considerably lower than 

A e ni| 7= 0- 

We then give a possible prediction for the entanglement 
for values considered reasonable by HAWH, by selecting 
the choice e = e opt given in their Fig. 2, and using the 
solution (|B3|> . The result is plotted in Fig. [3K and shows 
the possibility to detect entanglement in the presence of 
thermal mechanical decoherence, without the need to use 
laser cooling to reduce the value of no. 

Steering is confirmed if we satisfy the sum criterion ([3]) 
which becomes 

E ste er,sum(M\C) = 2{A 2 (A^ l " t + g x P° ut ) 

+A 2 (P™ t +g p X° ut )}<l. 

(B5) 

Thus, we can write the left side variances once again 
in terms of the steering parameter for 7 = 0, and the 



15 



thermal noise term, and normalise to obtain the measure 
of steering in the presence of mechanical decoherence: 

E steer . sum (M\C) = E gtSum (M\C)\ y= o + (2n + l)e. 

(B6) 

The entanglement criterion is symmetric with respect to 
the two systems, whereas the steering one is not. There- 
fore, we also calculate 

Al f (X° ut + g x P™*) = A 2 (X C + g x P™% =0 

+ \gl<n+\), (B7) 

^n S (P° ut + 9PKT) = a 2 (p c + 9 P X™%= (B8) 

+ \gl<n+^). (B9) 

Steering (C\M) is confirmed if we satisfy the criterion §3§ 
which becomes 

E steer , aum {C\M) = 2{A 2 (X° ut +g x P™ t ) 

+A 2 (Pr t + g P X° m ut )}<l, 

(BIO) 



and then, in this case, 



E st eer,8um(C\M) = E g , sum (C\M)\ 1=0 + {n + ^Hg 2 x + g 2 p ), 



where the g x and g p are optimised again, to minimise the 
variances and the measure of E ste er,sum(C\M). 

We plot the optimal steering prediction in Fig. [5J for 
the parameters chosen by HWAH in their Fig. 2. Of 
course, the observation of A ent < 0.5 is sufficient to verify 
steering, and this prediction can be read directly from 
the Fig. 2 plots presented by HWAH. However, A ent < 
0.5 will show steering only for g x = g p = 1, and we 
present the actual prediction for the steering parameter, 
optimised for g in Fig. 03. As for entanglement, the 
results indicate the potential to detect steering in the 
presence of thermal mechanical decoherence, without the 
need to use laser cooling to reduce the value of no- 
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